Mon. Not. R. Astron. Soc. 0Q0.mfl4l(2Q10) Printed January 6, 201 1 (MN KTeX style file v2.2) 



Lucky Imaging Adaptive Optics of the brown dwarf binary 
GJ569Bab* 

B. Femema^'^t, R. Rebolo^'^ J. A. Perez-Prieto^ S. R. Hildebrandt^''^, L. Labadie^'^, 
O A. Perez-Garrido^, V. J. S. Bejar^'^, A. Diaz- Sanchez^, 1. Villo^, A. Oscoz\ R. L6pez\ 
L. F. Rodriguez^ and J. Piqueras^'^ 



^Instituto de Astroflsica de Canarias, Q Via Ldctea S/N, E-38200 La Laguna, Spain 

^ Departamento de Astroflsica, Universidad de La Laguna, E-38205 La Laguna, Tenerife, Spain 
if 



Consejo Superior de Investigaciones Cientificas, Spain 
' Laboratoire de Physique Subatomique et Cosmologie, F-38000 Grenoble, France 

' ^ Universidad Politecnica de Cartagena, Campus Muralla del Mar, Cartagena, Murcia E-30202, Spain 

^ Max-Planck-Institut flir Sonnensystemforschung, Max-Planck-Str 2, D-37191 Katlenburg-Lindau, Germany 



Accepted 2010 December 16. Received 2010 December 15; in original form 2010 August 2 

Oh 



O 



(N 
> 

(N 
O 



X 
J3 



ABSTRACT 

The potential of combining Adaptive Optics (AO) and Lucky Imaging (LI) to achieve high 
precision astrometry and differential photometry in the optical is investigated by conducting 
observations of the close 0"! brown dwarf binary GJ569Bab. We took 50000 /-band images 
with our LI instrument FastCam attached to NAOMI, the 4.2-m William Herschel Telescope 
(WHT) AO facility. In order to extract the most of the astrometry and photometry of the 
GJ569Bab system we have resorted to a PSF fitting technique using the primary star GJ569A 
as a suitable PSF reference which exhibits an /-band magnitude of 7.78 + 0.03. The AO-hLI 
observations at WHT were able to resolve the binary system GJ569Bab located at 4"92 + 
0''05 from GJ569A. We measure a separation of 98.4 ± LI mas and /-band magnitudes of 
13.86 ±0.03 and 14.48 + 0.03 and / -/colors of 2.72+0.08 and 2.83+0.08 for the Ba and 
Bb components, respectively. Our study rules out the presence of any other companion to 
GJ569A down to magnitude I~ 17 at distances larger than 1". The I - J colors measured 
are consistent with M8.5-M9 spectral types for the Ba and Bb components. The available 
dynamical, photometric and spectroscopic data are consistent with a binary system with Ba 
being slightly (10-20%) more massive than Bb. We obtain new orbital parameters which are 
in good agreement with those in the literature. 

Key words: Instrumentation: high angular resolution - Instrumentation: adaptive optics - 
Stars: low-mass, brown dwarfs - Binaries: close - Stars:individual: GJ569 



1 INTRODUCTION 

The lucky imaging (LI) technique proposed by iFried ( Il978h at- 
tracted the attention by professional astrono mers once low read- 
out noise detectors became avail able (e.g. [Baldwin et alj [2OO1I : 
iTubbs et al]|2002l ; iLaw et alj|2006l) . Recently it has been realized 
that the combination of LI and Adaptive Optics (AO) can benefit 
mutually and provide high-resolution imaging close to the diffrac- 
tion limit at optical and near-infrared (NIR) wavelengths (e.g. see 

* Based on service observations made with the WHT telescope operated 
on the island of La Palma by the Isaac Newton Group and on observations 
made with the Nordic Optical Telescope, operated on the island of La Palma 
jointly by Denmark, Finland, Iceland, Norway, and Sweden, in the Spanish 
Observatorio del Roque de los Muchachos of the Instituto de Astroflsica de 
Canarias. 

t E-mail:bfemenia@iac.es 



lOladvsz et al.l2008l : lLaw et al.l2009l : [Kervella et alj2009l) . This can 
be obtained by using a larger fraction of data in the LI selection 
and/or keeping images with a better Strehl ratio than in conven- 
tional LI observations. 

In this paper we present the results of combining the LI 
and AO techniques to produce high-angular resolution and high- 
contrast imaging in the optical of the multiple system GJ569Bab 
which is a benchmark in the study and characterization of brown 
dwarfs (BDs). GJ569A is an M2.5V cr omospherically active 
star lying at a d i stance of 9.6-9.8 p c jPerrvman et al.l 1 19971 : 
Ivan Leeuwenll2007l) . lForrest etall fT98 s) identified a faint compan- 
ion to GI569 and argued the potential bro wn dwarf nature of s uch 
companion. Using Keck AO observations iMartm et al.l j2000l) re- 
solved GI569B as a binary brown dwarf system with a separation of 
~ 0" 1, a total mass of the system in the range 0.09-0. 15 M0 , an age 
in the range 0.12-1.0 Gyr and an orbital period ~ 3 yr. Further AO- 
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based observations with the Keck telescope (e.g . Lane et a0200ll: 
Zapatero Osorio et alj |2004| ; Isimon et alj l2006l: Konopackv et al 
2010) and the Subaru and HST telescopes fe apatero Osorio et al 



_2005; Martin et al. 2006) have allowed precise determination of 
the dynamical masses and orbital parameters of the binary system 
GJ569Bab (see Table[7]in this paper) as well as a precise determi- 
nation of the spectral types of the GJ569B component s: M8.5-9V 
and M9V for the Ba and Bb components, respectively jLane et al] 
l200lUMartin et al.ll2006h . 

For BDs an estimate of the mass is essential to determine their 
properties and evolution. A way to achieve a direct measurement 
of masses is to observe close binary systems where the short or- 
bital period allows for a complete sampling of the orbit and from 
here a precise determination of the dynamical masses of the pair. 
Up to now a complete characterization of the orbital motion has 
been ach ieved for a few BD binaries (e.g. see IDupuv et al.ll2009l : 
[Konopac kvetai] |2010|) . In this context GJ569Bab constitutes a 
unique laboratory where to test the stellar evolutionary models as 
it is among the shortest known period BD binary system. This has 
allowed to determine the orbit of the system over several periods 
and from here a precise determination of its dynamical mass. 

The high angular resolution requested to spatially resolve faint 
systems like GJ569B into its components has been so far achieved 
with 10-m class telescopes and AO in the NIR or with the HST. 
Our motivation to perform observations of GJ569 with LI-l-AO on 
a 4-m class telescope was twofold. First, to test the potential of this 
technique for high-angular resolution and high-contrast imaging in 
the optical regime. Second, to shed light on whether GJ569B is 
actually a triple syst em as suggested in some works in the liter ature 
jMartm e t al. 2000; Kenworthv et al.ll200ll : [simon et alj20o3) . 

In Section[2]we briefly describe the instrumental set-up. Sec- 
tion[3]describes the observations. Section|4]reports the data calibra- 
tion, reduction and analysis. Section [5] focuses on the photometry 
and astrometry of the GJ569 system components and discussion of 
the results. We provide our conclusions in Sect.|6l 



2 INSTRUMENTAL SETUP. 

FastCam is a Ll-based instrument developed jointly by the Insti- 
tuto de Astrofisica de Canarias (lAC) and Universidad Politecnica 
de Cartagena (UPCT). In brief, FastCam consists of a commer- 
cial Electro Multiplying Charge Coupled Device (EMCCD), a very 
simple versatile optics setup offering different plate scales to be ac- 
commodated at different telescopes configurations and an FPGA- 
basec{3 on-line processing and acquisition hardware. Its high de- 
gree of versatility has allowed us to collect images on a wide range 
of telescopes with diameters from 1.5 to 4.2 m. Further details of 
the FastCam instrumen t and results achieved so far at diffe r ent tel e- 
scopes are presented in lOscoz et al.l j2008l) ; lLabadie et"ar 

The observations reported in this paper were acquired at the 
4.2 meter William Herschel Telescope (WHT) at Observatorio del 
Roque de Los Muchachos (ORM, La Palma). FastCam was in- 
stalled at NAOMI, the AO facility at WHT. It suffices here to say it 
is an AO system based on a 8 x 8 subaperture Shack-Hartmann sen- 
sor and the most prominent feature of the system is the use of a fully 
segmented deformable mirror. Details for NAOMI can be found in 
[Myers et al. (2003); Benn (2006). We also resorted to non-AO LI 
observations with FastCam at the Cassegrain focus of the Nordic 



Optical Telescope (NOT) at the ORM mainly for plate scale and 
orientation calibration. 

The standard data acquisition of FastCam is that raw data are 
composed of cubes of 1000 images each. The full detector of Fast- 
Cam is an E2V EMCCD with 512x512 pixels (commercial camera 
model Andor iXon DU-987) exhibits a maximum frame rate of 35 
Hz, although the exposure time for each individual frame can be set 
at values shorter than 1/35 seconds but at the cost of incurring on 
overheads. 



3 OBSERVATIONS. 

3.1 AO-hLI observations with the WHT telescope. 

On the 4th of June 2009 FastCam was configured for an individual 
frame exposure time of 30 msec and with optics providing a field 
of view for FastCam of 7x7 arcsec^ (see Section |4jT3}. The pri- 
mary star GJ569A was used as the guide star for both the NAOMI 
and the FastCam real-time selection algorithm. Such real-time se- 
lection is performed by the FPGA-based acquisition system and it 
is mostly used to detect any anomalies or malfunctioning of the 
system by displaying the best 10% frames based on the brightest 
pixel selection although the entire frame set is kept for off-line pro- 
cessing. With this FastCam configuration, we collected 50 cubes 
which amount to a total exposure time of 1500 seconds in the stan- 
dard I Johnson-Bessell band (872.2 nm central wavelength, 298 nm 
FWHM) and 824 nm central wavelength, 175 nm FWHM when the 
EMCCD QE curve is considered. Taking into account overheads, 
the total time spent on GJ569 was 2340 seconds. During most of 
the observation time on GJ569 we had access to the RoboDIMlvj^ 
seeing monitor which reported an average seeing in the visible of 
0"49±0f' 08 (minimum seeing 0"38, maximum seeing 0"77) dur- 
ing the observation of GJ569. We also observed astrometric bina- 
ries for plate scale and orientation calibration (see Section l4. 1.2t . 



3.2 LI observations with the NOT telescope. 

On the 24th July 2008 FastCam was configured for an individ- 
ual frame exposure time of 20 msec and with optics providing a 
field of view for FastCam of 16x16 arcsec^ (see Section B.l.lt . 
In the I Johnson-Bessell band (same filter as with the WHT, see 
Section lSTI l we collected 10 cubes which correspond to a total ex- 
posure time of 200 seconds. No independent measurements for the 
seeing conditions were available at the time although the combined 
image of all 10000 individual frames yields an estimate for the see- 
ing of ~ 0"5 at the detector plane. The primary star GJ569A was 
used as the guide star for the FastCam real-time selection algo- 
rithm. We also observed the M15 globular cluster for plate scale 
and orientation calibration (see Section B.l.lt . 



3.3 Observations with the IAC80 telescope. 

For photometric calibration purposes we obtained /-band images of 
GJ569 on April 6 2010 using the CAMELOT instrument mounted 
on the IAC80 telescope with an identical filter to the ones used 
in Sections [37T] and [521 CAMELOT consists of a 2Kx2K CCD 
detector with a 0.304"/pix plate scale providing a field of view of 
10.4x10.4 arcmin^. We observed a series of 10 images with indi- 
vidual exposure time of 5s. Raw data were reduced using routines 



' FPGA, Field Programmable Gate Arrays 



' |http : / /www ■ ing . iac. es/Astronomy/development/seeing/Correl_DmM . htiiil| 
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Figure 1. Left: I-band image of GJ569 observed on the 4th of June 2009 with FastCam at the science focus of the WHT telescope AO system. Right: GJ569 in 
the I-band observed on the 24th of July 2008 with the NOT telescope. The image at the WHT is clearly resolved into the GJ569Ba and GJ569Bb components, 
while the NOT image acquired nearly 1 year before suggests a binary system because of the elongation of GJ569B. Both images have been rotated by their 
respective Position Angle Offsets (see Section l4~l1 to have the North and East axis pointing upwards and to the left, respectively). 



within the IRAF environment 0. Images bias was subtracted us- 
ing the overscan region and zero exposure time images. The flat 
field correction made use of dome flats. Aperture photometry was 
performed using routines within the DAOPHOT package. Weather 
conditions during our observations were photometric and average 
seeing ranged from 1.3 to 2". Instrumental magnitudes were trans- 
formed into apparent magnitudes in the Johnson-B essel system us- 
ing photometric standard stars from lLandol^(ll992h obtained at dif- 
ferent airmasses before and after our science target. Since images 
of GJ569 were defocused to avoid saturation, aperture correction 
was applied to the photometry of the standards to match the larger 
aperture used for our object. From these data, we have measured 
that the /-band magnitude of GJ569A is 7.78±0.03. 



4 DATA REDUCTION AND ANALYSIS 

The ima ge selection employed is the same as w ith non-AO obser- 
vations I Oscoz et alj|2008l : iLabadie et alj|20ld) which selects the 



best images as those exhibiting the largest ratio between the bright- 
est pixel in the image and the rms of the overall image. The images 
are sorted out from best to worst image according to the bright- 
est pixel criterion with a very small fraction of individual frames 
exhibiting saturated pixels (due either to cosmic rays or spurious 
electronic events) removed. Then the user defines the percentage of 



^ IRAF is distributed by National Optical Astronomy Observatories, which 
is operated by the Association of Universities for Research in Astronomy, 
Inc., under contract to the National Science Foundation. 



best images to be kept on which a shift-and-add algorithm is ap- 
plied. Standard image reduction considered only bias subtraction 
as experience with the FastCam EMCCD has revealed its excellent 
cosmetic and no need for flat fielding. Bias subs-traction, image se- 
lection and shift-and-add algorithm are performed using our own 
software. 



4.1 Plate Scale and Orientation Calibration 

4.1.1 NOT calibration with Ml 5 

As part of the observation program with the NOT telescope on the 
night of 24th July 2008, we also conducted observations of the M15 
globular cluster which, besides its specific science case, constitutes 
a very precise astrometric calibrator yielding a plate scale accu- 
racy of a few tens of microarcseconds and plate orientation below 
0.° 1 when the FastCam image of the M15 core is compare d against 
HST WFPC2 images of M15 jvan der Marel et alj|2002l) . For the 
instrumental setup on the 24th of July 2008 we measured a plate 
scale of 3 1.176 ±0.030 mas and a Position Angle (P.A.) offset of 
-89.84° ± 0.05° (i.e. North axis is at 89.84°clockwise from CCD 
y-axis). 



4.1.2 WHT calibration with WDS binaries. 

Observations of the M15 globular cluster were not possible on 
June 2009 as none of the cluster stars is bright enough for ef- 
ficiently locking the NAOMI loop. Because of this we observed 
four astrometric binaries (HD98231, HD105824, HD186858, 
HD197913) selected from the Washington Double Star (WDS) 
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Catalog jWvcoff etal.ll2006l) with an adequate component sepa- 
ration and brightness to act as astrometric calibrators. The cali- 
bration achieved from the WDS binaries (HD98231, HD105824, 
HD186858, HD197913) yields a plate scale of 13.63 +0.23 mas 
and a Position Angle Offset of 91.3° +0.9°. The uncertainties as- 
signed from the WDS astrometry corresponds to the dispersion of 
the values from each of the WDS binaries. 



4. 1.3 WHT calibration with GJ569. 

We are able to compute the GJ569B photocentre position with re- 
spect to GJ569A with high accuracy at the NOT image on July 24th 
2008 (see TableO. The orbital period of GJ569B around GJ569A is 
known a priori to be very long. With the information in Table|5]we 
predict the coordinates of GJ569B with respect to GJ569A on June 
4th 2009 (WHT observation) implicitly dismissing GJ569A-B or- 
bital displacements between the images at the NOT and at the WHT 
nearly one year later. From here a plate scale of 13.40 + 0.14 mas 
and a Position Angle Offset of 91.8° ±0.6° is derived. A weighted 
average of these values with those obtained in Section |4. 1 .2l vields 
the calibration assumed in this work for the WHT data: 13.46 ±0.12 
mas and 9 1 .6° ± 0.5° for the plate scale and Position Angle, respec- 
tively. 



4.2 PSF of the observations. 

The NAOMI system is optimized to deliver well corrected wave- 
fronts above 1.0 /im, while in the optical the aim is not to deliver 
well corrected images (i.e. with a significant Strehl ratio) but sim- 
ply to achieve a reduction by a factor 2 of the FHWM with respect 
to the seeing-limited images. The combination of AO+LI allows to 
achieve FHWM reductions by at least a factor 5 with respect to the 
seeing-limited images as shown in Fig.|2]and Fig. [3] where we show 
the evolution of the FWHM and different cuts of the PSF versus the 
LI % selection, respectively. The PSF cuts shown in Fig. [3] corre- 
spond to the PSF as measured on GJ569A for both the NOT and 
the WHT images. 

The FWHM evolution depicted in Fig. |2] shows the usual low 
sensitivity of the FWHM upon LI %, a tr end already pointed out in 
previous works (e.g. see iLaw et alj|2009l) . To understand this trend 
one must take into account that most of the turbulence is concen- 
trated in the tip-tilt modes which are removed in the LI scheme, 
while higher-order turbulent modes tend to distort the PSF but 
without enlarging it significantly. The effect of these higher-order 
modes in the PSF can be understood as a leak of power towards 
outside the PSF core, giving rise to the familiar LI PSFs consist- 
ing on a very well defined narrow core on top of a much wider and 
swallower halo. Qualitatively one can think that the better the AO 
correction or the most stringent the LI frame selection are, the more 
prominent the core over the halo becomes. This is in fact what we 
observe in the PSF cuts in Fig. [3] where the PSFs have been nor- 
malized to have unit volume. In this way we notice that when we 
restrict ourselves to very low LI % selection the PSF peaks much 
higher with respect to the halo that when using a larger LI % se- 
lection. This leads directly to a trade-off when having to decide 
the frame % to keep when one also considers the desirable SNR to 
achieve. 



PSF FWHM fvolutic 
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Figure 2. Evolution of the PSF FWHM at NOT (red) and WHT(black) as 
a function of the percentage of images kept in the selection procedure. The 
PSF is directly the image of GJ569A. In this plot the x-axis corresponds to 
the direction at which the PSF is widest and the y-axis is perpendicular to 
the X-axis. Notice these x and y axes do not necessarily correspond to the 
CCD x and y axes. The asymmetry of the PSF for the WHT observations is 
believed to be caused by the absence of an atmospheric dispersion corrector. 



4.3 The PSF fitting technique 

In order to perform an estimate of the relative photometry and as- 
trometry between the primary component and each of the resolved 
components of GJ569B we have resorted to a Le venberg-Marquadt 
(LM) non-linear least-squares l iBevingtonlll992l) PSF fitting algo- 
rithnfl We assume the PSF measured at the GJ569A location and 
the PSF at GJ569B position do not differ significantly. This is jus- 
tified by the excellent seeing conditions experienced during the ob- 
servations and the fact that the isoplanatic angle at the ORM at 
820 nm is frequently l arger than the 5 arcsec separation betw een 
GJ569A and B (see e.g. lFuensaUda et al.|l2004lLaw et al.ll2006h . It 



has also been pointed out bv lLaw et al.l 1 120061) that as the LI tech 



nique selects the best fraction of frames where the turbulence is 
smallest, it is expected that those frames also exhibit larger iso- 
planatic angles than would be observed in traditional long expo- 
sure imaging. Therefore, for a separation of about 5"we should be 
within the isoplanatic angle at ~ 820 nm and the assumption of 
uniformity of the PSF at the GJ569A and GJ569B locations is well 
justified. 

Our fitting model is applied to a relatively small box with only 
GJ569B and assumes the signal on each pixel in the box, /(x,y), is 
due to two components GJ569Ba and G569Bb located at positions 
(xa,ya) and (x/,,yi), respectively, and with fluxes a and b relative to 
the flux of GJ569A, plus a possible offset (OFF) associated either 
to sky background or a bias offset in the final image: 

f(x,y) = a-PSF(x-Xa,y-ya) + b-PSF(x-Xb,y-yt) + OFF (1) 

Note that in our LM PSF fitting scheme we do not fit an an- 
alytical PSF function nor attempt first a fit to the PSF measured 
at GJ569A. Instead we use directly the measured GJ569A PSF in 
the WHT image which is moved across a 100 x 100 pixel box by 
performing FFT-based shifts. In this way, the PSF fitting technique 
also provides us with a measurement of the distance between each 



* Using the MPFIT toolset by C. Markwardt available 
[http : //cow .physics . wise . edu/~craigm/idl/idl . htffll| 
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Figure 3. Cuts of the PSFs at the NOT and WHT as a function of the percentage of images kept in the selection procedure. The PSF is directly the image 
of GJ569A. The PSFs have been normalized so that the highest PSF exhibits peak unity and all of them have the same volume under the PSF.The gain in 
resolution from the WHT with respect to the NOT is due both to the increasing WHT diameter plus the use of the NAOMI system. Notice the blue line in 
the NOT plot represents the seeing-limited PSF while in the case of the the blue PSF in the WHT diagram corresponds to a partially tip-tilt removed by the 
NAOMI system. 



of the GJ569B components with respect to each other and with re- 
spect to GJ569A. 

The results of the LM PSF fitting technique on the image com- 
posed by the best 16% of the frames from the WHT observation are 
summarized in Table [T] for different values of the threshold back- 
ground parameter (rj). This parameter determines the number of 
pixels within the 100 x 100-pixel box which are finally considered 
in the fitting process: only those pixels exceeding the background 
level by a factor rj the rms background. Pixels below such thresh- 
old are shown in black in Fig.|4] As it can be seen in Table [T] when 
applied to the WHT image the technique is rather insensitive to the 
choice of pixels in the fitting (i.e. little sensitivity to the r] parame- 
ter) and provides the relative photometry with excellent accuracy. 

For the PSF fitting model to work properly and provide with 
meaningful error bars, the pixels in the image have to be assigned 
realistic error bars. During the analysis of the data the simple choice 
whether a given pixel was background limited or source limited was 
not good enough. At this point it revealed particularly useful the 
fact that the LI scheme allows for a precise distribution of values of 
the signal for each pixel in the final image since we have access to 
all the individual frames. For the final WHT image made of the best 
16% individual frames, the value assigned to each pixel is the mean 
of a distribution sampled with 8000 values and we also generated 
an error map computed as the error of the mean. 

The LM PSF fitting technique was also robust versus the LI se- 
lection fraction, exhibiting little dispersion on the derived param- 
eters as it is shown in Table |2] for the choice 7/ = 3.0. The choice 
of the 16% fraction to generate the final LI image for the WHT 
data was based entirely on the value of the reduced achieved 



in this series of fits: the WHT image generated with a 16% of the 
best images exhibits ^^/y = 1.04 (v being the number of degrees of 
freedom). 

The case for the NOT image is less favourable. When applied 
to the NOT image the LM PSF fitting was not that robust and we 
believe this is due both to the lower number of images under consid- 
eration but mainly to the fact the blending of both sources hampers 
a good fit. 



4.4 Detectability curves. 

A key point to characterize the combination of AO-assisted LI sys- 
tems for the detection of low-mass companions is the study of the 
detectability curves as a function of the distance to the star. This 
has been done in Fig. |5] where on the left panel we show the 3cr 
detectability curves for a variety of percentages of the images kept 
to generate the final images. On the right panel we have plotted the 
detectability curves on the same images upon filtering them with 
a wavelet transform algorithm. The details of this wavelet-based 
post-processing algorithm are described in Labadie et al. ( 201(3) 
and the general purpose of such filtering technique is the increase in 
the contrast ratio achievable by suppressing the low-frequency halo 
of the pr imary star which hides th e presence of much fainter com- 
panions jMasciadri & Ragal2004h . Both sets of detectability curves 
were obtained as follows: at a given distance p from the primary 
star we identify all possible sets of small boxes of size the FWHM 
of the PSF (i.e. 7x7 pixel boxes). Only regions of the image show- 
ing structures easily recognizable as spikes due to diffraction of the 
telescope spider and/or artifacts on the read-out of the detector are 
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Table 1. LM PSF fitting results on the WHT image from best 16% of frames. 





axlO^ 




bja 


p(inas) 


p.A. n 




2.0 


3.743 ±0.023 


2. 126 ±0.022 


0.568 ±0.009 


98.6±1.1 


78.0 ±0.7 


0.92 


2.5 


3.733 ±0.024 


2. 113 ±0.023 


0.566 ±0.009 


98.6±1.1 


78.0 ±0.7 


0.95 


3.0 


3.715 ±0.024 


2.096 ±0.023 


0.564 ±0.009 


98.4±1.1 


78.0 ±0.7 


1.04 


3.5 


3.689 ± 0.024 


2.076 ±0.023 


0.563 ± 0.009 


98.2±1.1 


78.0 ±0.7 


1.10 


4.0 


3.666 ±0.025 


2.054 ±0.024 


0.560 ±0.009 


97.8±1.1 


78.1 ±0.7 


1.16 


5.0 


3.626 ±0.026 


1.994 ±0.025 


0.550 ±0.009 


97.5±1.1 


78.2 ±0.7 


1.27 


6.0 


3.561 ±0.028 


1.920 ±0.026 


0.539 ±0.010 


96.8±1.1 


78.1 ±0.7 


1.33 


7.0 


3.490 ±0.029 


1.821 ±0.028 


0.522 ±0.010 


95.8±1.2 


78.1 ±0.7 


1.25 



Table 2. LM PSF fitting results as a function of the percentage of frames kept to generate the WHT image. 



LI% 


axlO' 




h/a 


p(mas) 


P.A. n 




1 


3.35 ±0.06 


1.77 ±0.05 


0.528 ±0.021 


96.5 ±1.8 


77.8±1.1 


0.880 


5 


3.59 ±0.05 


2.07 ±0.04 


0.576 ±0.018 


98.1 ±1.5 


78.2 ±0.9 


0.360 


10 


3.688 ±0.029 


2.082 ±0.026 


0.565 ±0.010 


98.5 ±1.2 


78.2±0.7 


0.722 


16 


3.715 ±0.024 


2.096 ±0.023 


0.564 ±0.009 


98.4±1.1 


78.2 ±0.6 


1.04 


20 


3.7 10 ±0.022 


2.125 ±0.021 


0.573 ±0.008 


98.3±1.1 


78.0 ±0.6 


1.27 


40 


3.745 ±0.017 


2.148±0.016 


0.574 ±0.006 


98.3 ±1.0 


78.3±0.6 


2.11 


60 


3.765 ±0.015 


2.140±0.014 


0.569 ±0.005 


98.9 ±1.0 


78.2±0.6 


2.77 


80 


3.755 ±0.014 


2.132±0.013 


0.568 ±0.005 


98.6 ±1.0 


78.1 ±0.5 


3.24 



not considered. For each of the vahd boxes on the arc at distance p 
the standard deviation of the image pixels within the 7-pixel boxes 
is computed. The value assigned to the 3cr detectability curve at dis- 
tance p is 3 times the mean value from the standard deviations of 
all the eligible boxes at distance p from GJ569A. From the curves 
in both panels in Fig.[5]we notice the following behaviours: 

(i) On the non- wavelet filtered image, at short distances (< 1"5) 
from the primary star, the lower the fraction of data used the higher 
the detectability gain is with a maximum gain of around ~ 0.9 
magnitudes between the 80% and 1% frame selection. This im- 
provement in contrast is due a slightly better angular resolution 
with a very restrictive frame selection. As we move away from the 
primary star, the detectability curves approach to each other until 
eventually the situation is reversed. At large distances the image 
is dominated by background/detector noise so that the higher the 
percentage of frames used to generate the image, the deeper the 
detectability curve gets. 

(ii) On the wavelet-processed image there is no transition be- 
tween the inner region (i.e separation from primary below 1"5) 
and the outer region. The detectability at the inner regions is always 
improved with respect to the non-processed image and, in addition, 
the higher the percentage of data kept for the image, the deeper 
the detectability gain. Above If 5 separation from the primary star, 
there is no real difference between the detectability achieved on the 
original and the wavelet-processed image, as it is expected at those 
regions of the image which are dominated by background/detector 
noise which can not be removed by the wavelet post-processing 
technique. 



5 RESULTS. 

5.1 Photometry of the GJ569 system. 

The final /-band differential magnitudes (AniAB = ms - niA) de- 
rived by the LM PSF fitting technique described in Section 14.31 



Table 3. Differential magnitudes in I band with FastCam at WHT 



AifiAB 5.590 ±0.004 

AniABa 6.075 ±0.007 

AniABb 6.697 ±0.012 

^niBaBb 0.622 ±0.017 



are summarized in Table [3] The magnitude difference between 
GJ569Ba and GJ569Bb we find is very similar to the 0.7 mag 
found bv Marti n etai] ( 1200 6^ using the F814W filter in the HST 
which resembles the /-band. These values are slightly larger than 
the 0.5 magnitude difference measured in the J,H and A'-band by 
iMartm et al.l JioOQ) . All these data are consistent with the spectral 
determination of M 8.5-M9 for t he bro wn dwarf binary derived in 
the near-infrared bv Lane et al.l ( l200ll) . However, given the uncer- 
tainties, we can not discard that b oth components hav e a similar 
spectral type of M9, as obtained bv lMartm et al. ||2003) using opti- 
cal spectra. 

Using the /-band photometry obtained in the IAC80 for 
GJ569A (see Section [331 l, and the differential magnitude between 
the primary and the brown dwarf binary found with FastCam (see 
Table|3), we derived an integrated photometry of /=13.37±0.03 for 
GJ569B with the WHT image and /=13.41±0.03 with the NOT im- 
age. As remarked previously, the lack of enough spatial resolution 
on the NOT image does not allow to resolve the GJ569B binary 
system, but aperture photometry of the whole GJ569B system is 
compatible with what we derive from the LM PSF fitting of the 
WHT im age. This valu e is not in good agreement with the previous 
resuh bv lForrest etaT ( 1988) (/=13.88±0.2), probably affected by 
the contamination of the primary. To check a possible variability 
of GJ569A , we have compared the apparent magnitude of this ob- 
ject and the GIV star HD130948, which was observed during the 
same nights in the IAC80 and the NOT. We have found no relative 
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GJ569Bab with WHT on 04/06/2009 



Data vs. Best Fit 




Figure 4. Comparison of actual images to best fit models and respective cuts along lines joining GJ569Ba to GJ569Bb for data acquired with the WHT on 
4/6/2009. The oversampling of the WHT data and a better PSF due to the use of the AO facility at the WHT allows to resolve the binary nature of GJ569B 
allowing the LM PSF fitting technique to achieve a good fit to the data. On both images dark areas indicate regions whose pixels are not sufficiently above the 
background threshold to be taken into consideration in the PSF fitting (e.g. parameter in Tables [T] and [2). No rotation to aUgn the North and East axis with 
the usual y and negative axes has been applied to the image. 



3 — CT Detectability vs LI selection 




1 2 3 4 1 2 3 4 



Distance to GJ569A [arcsec] 

Figure 5. Evolution of the 3<t detectability versus the percentage of images kept to generate the LI image. The left panel shows the 3(T detectability curves 
obtained directly from the image. Right panel shows the 3-o" detectability curve upon filtering with wavelets (see main text). The green thick solid lines in 
both plots correspond to a selection of the best 16 % frames without wavelet filtering and is displayed in the right panel for comparison purposes. Notice the 
big improvement by the wavelet filtering at angular distances below 1"5. a slight gain at angular distances in the range [lf'5, 2"] and no net gain beyond 2". 
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Table 5. Astrometry of the binary system GJ569A-GJ569B. 



10 

Sp. Type 



15 



Figure 6. / — 7,spectral type diagram of GJ569Ba and GJ569b (solid red 
circles) in comparison with M (M5 to M9 coiTespond to x-axis values 5 
to 9) and L field dwarfs (L P to L8 correspond to x-axis values 10 to 18) 
from lLiebert & Gizij <200d) (open circles). The average color of these field 
dwarfs with the same spectral t ^pe are ind icated by the solid line. The J- 
ban photometry obtained from Lane et alj bOOlh in the CIT s ystem has 
been t ransformed into the 2MASS system using the relations in lCarpenteJ 
<200ll) . 



difference in their /-band magnitude at the level of a few tens of 
millimagnitudes. The absolute photometry of both components of 
the brown dwarf binary was conducted in a similar way and a sum- 
mary with the apparent magnitudes as derived from this work and 

in the literature is given in Table |4] 

With the 7-band photometry in lLane et alJ ( 1200 l l) we derived a 
1-J color of 2.72±0.08 and 2.83+0.08 for GJ569Ba and GJ569Bb, 
respectively. In Fig. [6] we show the / - / color of both components 
in co mparison with M and L field dwarfs of the same spectral type 
from lLiebert & Gizij ilOOdt) . We have adopted the spectral deter- 
mination of M 8.5-M9 for t he bro wn dwarf binary derived in the 
near-infrared bv lLaneet ail toll) . From this figure we can see that 
our estimated I- J colors are in good agreement with these spectral 
types. 

Ac cording to the absolute magnitudes and colors of field 
dwarfs ('Kirkpatrick & McCarthy' 1 994 iLiebert & Gizisll2006l) . we 
expect a relative magnitude difference of ~ 0.4 mag in the /-band, 
and a / - 7 color difi'erence of 0.27±0.22 for a M8.5/M9. All these 
data are also consistent with previous determinations of spectral 
types, in particular for those deriving a slightly earlier spectral type 
for GJ569Ba than for the Bb component. 



5.2 Astrometry of GJ569A-GJ569B. 

The compilation of all astrometric data on the GJ569B orbit around 
GJ569A, including our new measurement from the NOT observa- 
tion, is provided in Table|5]and shown in Fig.|7|with the value ob- 
tained in this work in red. The higher accuracy of our measurement 
is immediately noticeable and it is due to the very precise image 
statistics combined with the high angular resolution PSFs provided 







Separation 


PA 


Date 


Ref. 


(arcsec) 


n 


1985 Jul 28 


U) 


5.07±0.15 


17.4±1.7 


1987 Sep 18 


(a) 


5.01 ±0.15 


17.3 ±1.7 


1999 Aug 29 


(b) 


5.00 ±0.05 


25.0 ±0.9 


2001 Feb 8 


(c) 


4.89 ±0.04 


30±3 


2005 Feb 25 


(d) 


4.99±0.19 


33.1 ±1.9 


2008 Jul 24 


(e) 


4.94 ±0.04 


34.90 ±0.05 



Refere nc es: ^"jForres t et ^ 
i200ll) :"' ISh^on et alj feOOQ) : 



h988h : ' NMartin et al.l <2000l) : '' iLane et alJ 

'■''This work 



Oxbit of CxJ569B around GJ569A 



Sep 198? 



Aug 1999 







Feh 2001 










Feb 2005 






200B 





5.0 ^ 
4.8 - 

- 4.6- 

m 

o 

^ 4.4- 
o 

c 4.2^ 
4.0 

3.8 [ , , , \ , , , , \ , , , , \ ^ 

-2.5 -2.0 -1.5 

RA X cos(decl.) (arcsec) 

Figure 7. Compilation of positions of GJ569B on its long-period orbit 
around GJ569A using measurements in Table \5\ and our point in red. 
GJ569A would be located at coordinates (0,0), which falls on purpose out- 
side of the plot in order to make visible all points with their respective un- 
certainties. 

by the LI technique and a very accurate NOT plate calibration, as 

discussed in Section|4T| 

Only the photocentre positions in lForrest et al.l ( Il988l) and in 
this work correspond to observations in the I band, while the rest of 
the observations have been conducted in the NIR. Since the flux ra- 
tios between the GJ569B components differ slightly from the NIR 
to the I-band measurement, we realize there might be a slight bias 
in the photocentre determination. In principle such a bias could be 
avoided by a proper weighting with the individual masses of the 
GJ569B pair, but due to the uncertainties in the individual GJ569Ba 
and GI569Bb masses, the estimate of the barycenter remains essen- 
tially the same as the computation of the photocentre. The limited 
coverage of the GJ569B orbit around GJ569A (see Fig.|7]l does not 
allow for a reliable orbit computation but still we estimate a low 
eccentricity long-period orbit of ~400 years with a semi-major axis 
of 4f'3. 



5.3 Orbital analysis of the GJ569B system. 



As discussed in Section |431 the LM PSF fitting procedure applied 
to the WHT images allows us to obtain with high precision the rel- 
ative position of Bb with respect to Ba (see Table|6l(. The same PSF 
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Table 4. Updated compilation of GJ569 photometry. 





Spec. Type 


V magnitude 


I magnitude 


J magnitude 


H magnitude 


K magnitude 


GJ569A 


M2.5V<* 


10. lit*''' 


7.78±0.()3<J> 


6.633 ±0.023<'=> 


5.990 ±0.021<'=' 


5.770±0.018W 


GJ569B 
GJ569Ba 


M8.5-9V(''-f> 




13.88 ±0.2<"' 
13.37 ±0.03*® 
13.41 ±0.03"<J' 
13.86 ±0.03'J' 


10.61 ±0.05<* 
10.75 ±0.06(" 

11.14±0.07<* 
11.18±0.08W 


10.16±0.10<''> 
10.15±0.04<'> 

10.43 ±0.O4<!5> 


9.56±0.10'"> 
9.45 ±0.05'* 
9.62 ±0.03® 
10.02 ±0.08<* 
9.86 ±0. ID'S' 


GJ569Bb 






14.48 ±0.03'J' 


11.65 ±0.07** 
11.69±O.08W 


11.04±0.05<s> 


10.43 ±0.08<* 
10.39 ±0.06<s> 


Bb/Ba 






0.70±0.231:<'' 
0.622 ±0.017tj' 


0.5±0.2<'^' 
0.51±0.02<* 


0.5 ±0.1^ 
0.61±0.03<s) 
0.57 ± 0.04^) 


0.5 ±0.1"^' 
0.41±0.03<* 
0.41 ±0.13® 
0.61±0.03's' 
0.500 ±0.008*' 



R eferences: '^Forrest et all jl988l); <''lPerrvman et all jl997l); ''^iMartmetal] j2000l) ; <'' lLane et alj i200ll) ; '" ICutri et alj )2003l) :' lM^in et al.l j2006l) ; 
te lSimon et all i2006h : "' iKonopackv et alj bOlOl) : '■ IDuduv et al.l j2010h :^'This work. 

Notes: (t)lO.l 1 V Johnson magnitude converted from Hipparcos photometric system Hp = 10.201 ±0.004; '**Unresolved GJ569B /-band magnitude in June 
2009 at WHT; '"'Unresolved GJ569B /-band magnitude in July 2008 at NOT; (t)Differential magnitude using the F814W filter in the HST which resembles 
the /-band. 



fitting approach applied to the unresolved GJ569B NOT image de- 
livers a relative astrometry for these two components with higher 
error bars. Fig. [8] summarizes all the observations in the literature 
with our measurements in red. Our WHT astrometric value on June 
20 09 is nearly over im pos ed on the Keck AO-base d observations 
bv lDupuv etai] ( l2010h and IKonopackv etalj feOlOl) . the latter ob- 
servved GJ569B within seven days after our observation. Overplot- 
ted we show the different orbits found in the literature together with 
our best Keplerian orbit to the same data plus our two new points. 
The orbital parameters of all these orbits are provided in Table|7]to- 
gether with the total mass of the GJ569B system obtained by direct 
application of Kepler's third law. A very good agreement is found 
for the orbital parameters as well as for the estimated total mass 
of the system except from the value derived bv lDupuv et alJdlOlOl) 
(see discussion later in this section.) 

In Table |7] in addition to the orbital parameters of the brown 
dwarf system, we also list the total mass estimates found in the lit- 
erature together with our estimation which is consistent but slightly 
lower than previous results. 

During the referee process of this work IDupuv et alj ilOld) 
published a new set of values for the GJ569B system astrom- 
etr y. As the aut h ors re mark, there has been a recent revision 
by Ivan Leeuwenl | |2007|) of the GJ569 Hipparcos parallax yield- 
ing a revised distance of 9.65 ±0.16 pc instead of previous es- 
timates of 9.81 ±0.16 pc jPerrvman et al.l [l997l) ; the impact of 
this new estimate of the GJ569 system distance is shown as 
an additional entry in Table I?] iDupuv et al. (2010) claim pre- 
vi ous astrom e tric v al ues obtained with th e Keck AO facility 
in iLane et alj ( 1200 ll) ; IZapatero Osorio et alj j2004l) ; ISimon et all 
( l2006l) assumed a wrong plate scale. The same situation happens 
fo r the HST value o n 2002 June 26th (which is also considered 
in lZapatero OsoritTet a l. 2005). Dup uy et al. (2010) reanalyze both 
the astrometric data in lsimon et al. I ( 12006 ) and the HST value on 
2002 June 26th, and together with new 5 data points from January 
2008 until May 2010 a new orbit is computed providing a mass 
estimate for the GJ569B system of 0.140+|] |][]^ M© . The com- 



Table 6. Astrometry of the binary system GJ569Ba-Bb. 







Teles- 




Separation 


PA 


Date 


Ref 


cope 


In fit 


(arcsec) 


(°) 


1999 Aug 29 


(a) 


Keck 


Yes 


0.101 ±0.001 


46.8 ±3 


2000 Feb 18 


(b) 


Keck 


Yes 


0.092 ±0.001 


98.2±3 


2000 Feb 25 


(b) 


Keck 


Yes 


0.090 ±0.001 


100.4 ±2 


2000 Jun 20 


(b) 


Keck 


Yes 


0.076 ±0.003 


138.6±2 


2000 Jul 4 


(c) 


MMT 


No 


0.078 ±0.003 


148 ±3 


2001 Mar 9 


(c) 


SAO 


No 


0.0896 ±0.0010 


321 ±1 


2001 Mar 10 


(c) 


SAO 


No 


0.0899 ±0.0010 


320 ±1 


2000 Aug 24 


(b) 


Keck 


Yes 


0.059 ±0.001 


178.4±2 


2001 Jan 9 


(b) 


Keck 


Yes 


0.073 ±0.002 


291.4±2 


2001 May 10 


(b) 


Keck 


Yes 


0.097 ±0.001 


341. 1±3 


2001 Jun 28 


(tl) 


Keck 


Yes 


0.1024 ±0.0012 


352.6±2 


2001 Sep 1 


(tl) 


Keck 


Yes 


0.1033 ±0.0010 


9.7 ±2 


2002 Jun 26 


(e) 


HST 


No 


0.090 ±0.003 


94.4 ±1.2 


2003 Jul 1 1 


(e) 


Subaru 


No 


0.085 ±0.007 


319.8±4.4 


2004 Dec 24 


(f) 


Keck 


Yes 


0.0885 ±0.0048 


111.3±1.2 


2005 Feb 25 


(f) 


Keck 


Yes 


0.0798 ±0.0029 


133.7 ±0.7 


2008 Jan 16 


(E) 


Keck 


No 


0.061 8 ±0.0007 


272.5 ±1.4 


2008 Jul 24 


(h) 


NOT 


Yes 


0.0981 ±0.0022 


351.8±1.4 


2009 Apr 29 


(E) 


Keck 


No 


0.1002 ±0.0006 


66.5 ±0.6 


2008 May 29 


(E) 


Keck 


No 


0.1005 ±0.0018 


75.0±1.7 


2009 Jun 4 


(b) 


WHT 


Yes 


0.0984 ±0.0011 


78.2±0.6 


2009 Jun 11 


(i) 


Keck 


Yes 


0.099 ±0.002 


79.4 ±0.2 


2010 Mar 22 


(E) 


Keck 


No 


0.0558 ±0.0003 


206.5 ±1.1 


2010 May 23 


(E) 


Keck 


No 


0.0599 ±0.0004 


268.1 ±0.4 



R eferences: '^iMiu^m et al.| ( 2000);"'- L ane et alj 

'■^ Kenworthv et alj I20qil) ;'''-Zapatero Osorio et alj 

; Zapatero Osorio et alj j2005l); '"Is imon et alj J2006l) ; 



2001): 



2004); 



IDupuv et alj 



bOlOl) ; ''''This work; " iKonopackv et al., qOlQ) 
Column "In Fit" flags the data points used in our orbital solution. Data 
points flagged with NO correspond to new additions appeared during the 
referee stage. 
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Table 7. Orbital Parameters of GJ569Bab. 



Parameter Zapatero04*'' Simon06''* KonopackylO'^' DupuylO*^' This Paper 



Semimajor axis, a (mas) 


92 ± 2<"> 


90.4 ±0.7 


90.8 ±0.8 


95 6+'-' 


90.1 ±0.7 


Eccentricity, e 


0.32±0.01 


0.3 12 ±0.007 


0.310±0.006 


0.316 ±0.005 


0.317±0.010 


Inclination, i (°) 


34±2 


32.4±1.3 


33.6±1.3 


35.0±1.1 


30.0 ±1.6 


Periapsis argument, oj(°} 


257 ±2 


256.7 ±1.7 


77.4±1.7<'" 


257.9 ±2.0 


262.5 ±2.2 


Ascending Node Longitude, 


321.5±2.0 


321.3±2.2 


144.8 ±1.9<'" 


324.8 ±2.0 


313±4 


Peiiod, P(days) 


876 ±9 


863.7 ±4.2 


865.1 ±0.7 


864.0 ±1.1 


870 ±9 


Total Mass, M(Mo 


0.125 ±0.005 


0.125 ±0.007 


0.1 26 ±0.007 




0.122 ±0.007 


Total Mass, M(Mo 


0.1 19 ±0.005 


0.1 19 ±0.007 


0.120 ±0.007 


140+"™' 


0.1 16 ±0.007 


Epoch, TCyrs)'^* 


2000.758 ±0.008 


2000.754 ±0.007 


2003. 150 ±0.005 


2010.255 ±0.004 




Epoch, TCyrs)'*^ 


2010.36±0.10 


2010.22 ±0.05 


2010.260 ±0.008 


2010.255 ±0.004 


2010.263 ±0.012 



References: " IZapatero Osorio et al.l i2004 . '^ ISimon et all i2006h . '^ iKonopackv et al.l i2O10l) , '"^ IDupuv et all j2010l) 

Notes: 

In original paper quoted in AU. Converted to angular distance assuming GI569 Hipparcos parallax n = 0.10191 ±0.00167 mas iPenTman et al J 19971) . 
oj and n d eiived solely from astrometry results exhib it a 1 80° degeneracy which is broken when using spectroscopic data (e.g. radial velocities) such as 
in lZapatero Os orio et al. (2004), Simon et al. (200^ and lKon opackv et al. (2010). However in the latter work there seem s to appear again the degeneracy 

which was broken in iZapatero Osorio et ah Ii2004l) . The reason why such a degeneracy remains in Konopacky et al. 1 2010) is unknown to us^ 

To t al Mass deriv ed using th e GJ569A Hipparc os parallax tt = 0.10191 ± 0.00167 mas (j^erryman et al. 1997) in original works bv .Zapatero Osorio et alj 
j2004l) . ISimon et alj 12006) and Konopackv et al.l ilOlO). 

'■'^ Total Mass computed with the revised Hi pparco s parallax ;r = 0. 10359 ± 0.00172 mas ivaii LeeuwerjIioOTi) as pointed out in IDupuv et alj j2010l) . The 
total mass estimates in lZapatero Osorio et aljj2004h . lSimon et alj j2006h and lKonopackv et^n20Kl) scaled to the new parallax value. 
Epoch of periapsis appearing in original works. 

Epoch of periapsis propagated to fall within 2010 for the ease of comparison. 



parison of the new orbit by ^ Dupuv et alj ( l20ld ) against previous 
orbits in the literature is shown in a fraction of the orbit span- 
ning the period April - June 2009 in Fig.|9] From theoretical mod- 
els £uguZstal] bOlQ) estimate the mass fraction parameter to be 
q ^ 0.87 ± 0.02, which agree v ery w ell with observational r esults 
in lZapatero Osorio et al J j2004l) and ("Konopacky et alj |2010l) . The 
new astrometric values in lDupuv et aL C2010.) have not been con- 
sidered in our orbital solution. 



5.4 Constraining tlie binary nature of GJ569Ba. 

From radial velocities IZapatero Osorio et al. I j2004l) assign m- 
dividual masses of 0.071±0.011 Mq and 0.054±0.011 Mq 
for the GJ569Ba and GJ569Bb components, respectively. This 
is consistent wi t h 0.() 73±0.008 Mq and 0.053±0.006 Mq by 
iKonopackv et alj ( 1201 Oh from radial velocity measurements. How- 
ever, upon the dis covery of the multiple nature of GJ569B by 
iMarti'n et ^ i2000l) . the possi bility of GJ569Ba b inarity was dis- 
cussed. The reason argued in i Martm et alj (12000*) was the failure 
of theoretical models to fit the observational data with a single 
isochrone. The b i nary n ature of GJ569Ba was again suggested in 
iKenworthv et alj j200lh on the basis of Bb being about half as 
bright as Ba while the GJ569B spectrum was fit ted by a singl e MB. 5 
component and without evidence of reddening. ISimon et alj | |200^ 
gives further support to the Ba binarity when accounting from their 
estimated mass ratio q =MBb /Mbb =0. 19±0.13 as derived from 
their radial velocities. ISimon et al. |(l200i) suggest GJ569B is actu- 
ally a triple system with Baa, Bab (i.e. the two suggested compo- 
nents of Ba) and Bb with nearly the same masses of around 0.04- 
0.05 M0 . H ow does such a pos sibility fit into the larger collection 
of data since lSimon et alj ( |2006|) ? 

(i) The flux ratio evolution with wavelength is not consistent 
with a system of nearly three-equal mass objects. In such a case the 
flux ratio should remain constant while its evolution from I-band 
to L-band is clearly decreasing: 0.622 + 0.017, 0.51+0.02, 0.590 + 



0.023, 0.504 + 0.015 ar id 0.49 ± 0.03 in I . J, H, K and L' bands 
(the L' measurement in iDupuvetalJboiOl) . This flux evolution is, 
however, consistent with Ba and Bb exhibiting slightly different 
radii (Ba larger than Bb) as suggested from theoretical models of 
BDs for the assumed masses of Ba and Bb and a slightly change 
in the effective temperature (Ba slightly hotter t han Bb) consisten t 
with the assigned spectral types to Ba and Bb in lLane et alj ( l200lh . 
Notice also that the a bsence of a Li feature in th e combined spec- 
trum of GJ569B in tZapatero Osorio et alj|2003) plays against the 
three-equal mass objects hypothesis: assuming the highest mass 
estimate for the GJ569B system of 0.140+° °°^ Mq dDupuv et al.l 
2010) would inevitably ensure the presence o f Li features as all ob- 
jects would be below the 0.060 M© t hreshold jChabrie r et aL 2000|) 
but not such a feature was detected in lZapatero Osorio et aL (20050 . 

(ii) The a ^e of the GJ569B is still an open issue with 
ISimon et al. I (2006^ assigning a similar age to the Pleiades (~ 
120 Myr) and the rest of estimates in the literature proposing 
ages in the range 200-800 Myr. The (1-J) colors of 2.72+0.08 and 
2.83+0.08 for the Ba and Bb components, respectively, do not 
match with what is expected from Pleiades objects with masses 
~ 0.04 - 0.05 Mq in the J v s. (1-J) Pleiades color-magnitude di- 
agram in iBihain et alj ( I2OO6I) . Note that the older the system the 
more inconsistent the measured (1-J) colors of Ba and Bb with re- 
spect to masses as lo w as 0.04-0.05 M p . Note also that according 
to the the Lyon (Cha brier et al. II2OOOI) and Tucson dBurrows et alj 
Il997h evolutionary models and for ages older than 250 Myr, ob- 
jects with masses in the range 0.04-0.05 Mq should exhibit spec- 
tral types in the L and T domain which is in contradiction with the 
M8. 5-9.0 spectral types assigned to the unresolved GJ569B and re- 
solved Ba and Bb components. 

However, there is still the possibility that Ba is composed of 
two components, namely Baa and Bab, of different masses. For 
the purpose of the following discussion on the restrictions on the 
GJ569Ba binarity we assume the tot al mass estimate of Mjot = 
0.140+[j{j[j^ Mq in lDupuv et si\ ilOVj) and a mass ratio q = 0.80 + 
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Figure 8. Compilation of astrometric data on the GJ569Bab pair witli our points in red (see Table [6). Overplotted there are the different best-fit Keplerian 
orbits for the GJ569 multiple system found in the literature and derived in this work and whose orbital parameters are summarized in TablefT] 



0.15 which h es midway and consis t ently betw een the observed 

! ' valij es in in lZapatero Osorio et alj ||2004) an d lKonopackv et alj 
2010 l). and the estimates bv lPupuv et aljfeoioh of 0.866+° °[^ and 



0.886 



+0.021 



using the Lyon and Tucson evolutinary models. With 



these Mxot and q values we have = 0.078 + 0.008 M© and and 
Mflb = 0.062 ± 0.008 Mq . A c lose look at the ast rometric values 
re ported in April-May 200 9 bv lPuDuv et al] feOlOh . by June 2009 
in'Konopackv et all j2010 f) and in this work are shown in Fig. [9l 
allowing us to place a constrain to the combination of mass and 
orbit semi-axis of a possible component Bab. From the above dis- 
cussion we may safely assume that if GJ569Ba is binary then its 



most massive component must satisfy 0.060 <MBaa [Mq < 0.086 
and <MBab /Mq < 0.018. Also, if Bab exists it should be ei- 
ther very faint and orbit around Baa with a semi-axis smaller than 
0.5 AU (i.e. midpoint between Ba and Bb) and, using Mga = 
0.078 ± 0.008 M0 , a possible orbital period T< 1 .3 yr. The presence 
of Bab would induce a motion in the position of GJ569Baa which 
would reflect itself as a wobble of the orbit of GJ569Bb around the 
compu ted Keplerian orbits and with an amplitude a given by (e.g. 
IPerrvm an 2000): 



Mp a 
M^d 



(2) 
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Figure 9. Zoomed part of orbit showing orbit in period April- June 2009 
to gether with orbits in the literature and predicted positions in orbit model 
bv lDupuv et"ZI iaOlOll . 



where Mp and are the masses of GJ569Bab and GJ569Baa, 
respectively, a is the Baa-Bab orbit semi-axis and d the distance 
to GJ569 (i.e. d ~ 10 pc). Examining in detail th e devi ation of ob- 
served positions with respect to the Dup uv et alj ( 1201 (j) orbit in the 
two month period depicted in Fig. [9] we place an upper limit a < 
2.2 mas with Icr confidence level. Using Msaa < 0.086 Mq we ob- 
tain the constrain a-MBab <1.93, with Mfiab in units of Mj (Jupiter 
mass) and a in AU. On the other hand we have placed an upper 
limit T< 1.3 yr on the period of Bab around of Baa. The resid- 
ual rms of the observed positions on the period April 2009-May 
2010 with respect the predicted positions allows to place an upper 
limit ff < 1.9 mas with Itr confidence level and from here the con- 
strain a-Msab <1.68. Adopting the former and more conservative 
constrain and since the minimum separation between GJ569Ba and 
GJ569Bb is ~ 0.5 AU, the widest possible orbit of Bab around Baa 
would be < 0.25 AU and we would have Meab ^ 8 Mj and the 
closer Bab to Baa the higher the feasible value for Meab • Given 
the range of masses under consideration and the distances between 
Baa and Bab, radial velocities studies of GJ569Ba should be able 
to shed further light into the existence of the Bab component. 



6 CONCLUSIONS. 

We have presented results showing the potential for high precision 
astrometry, differential photometry and high contrast imaging us- 
ing a Lucky Imaging instrument coupled to an adaptive-optics sys- 
tem. Our results indicate that with 4-m class telescopes equipped 
with a moderately low order adaptive optics system it is possible to 
achieve angular resolutions better than 0" 1 in the I-band. This is 



comparable to what is achieved in the Ks band with the use of AO- 
systems at 8-10 m class telescopes. This work is part of an effort 
to determine the feasibility of a Lucky Imaging instrument to be 
coupled on the future AO system at the GTC telescope (GTCAO) 
and what the expectations in terms of high contrast imaging and 
angular resolution should be expected from such a combination. 

Our work has focused on the observation of the GJ569 system 
which contains a benchmark brown dwarf binary. The application 
of our LM PSF fitting technique to the GJ569 image obtained with 
the WHT, where the GJ569B components are resolved, allows to 
achieve high precision relative photometry (to a few millimagni- 
tudes) and astrometry (to a few mas) thanks to the availability of 
meaningful error bars associated to each pixel in the fitted image. 

The potential for detection of faint companions has been ad- 
dressed by looking at the 3cr detectability curves in Fig.|5] On the 
images directly from the frame selection procedure, we distinguish 
two regions in which the detectability behaves differently versus 
increasing the percentage of images being employed. In regions 
where the image is dominated by the wide swallow halo of the 
primary PSF, the detectability is improved by restricting the per- 
centage of images being employed since in that way there is more 
energy in the core and less in the halo. At large distances from 
the primary, the image is dominated by background/detector noise 
and the detectability is improved by simply adding as many im- 
ages as possible. On the wavelet-processed images, both in the in- 
ner and outer regions, we only see a benefit on the increase of the 
percentage of images used in the frame selection. With the wavelet- 
processed images we observe a magnitude gain in the inner region 
(1.7 magnitudes at I"with respect to the non-processed image) but 
far away from the primary no net gain as those parts of the image 
are dominated by background/detector noise. 

We have measured a differential magnitude at I band between 
GJ569Ba and GJ569Bb hmgaBb = 0.622 ± 0.017. When used in 
conjunction of AniBaBb in the J, H and K bands by previous works, 
fits well with the spectral determination of M 8.5-M9 for the b rown 
dwarf binary derived in the near-infrared bv lLane et alj feoOlh. Our 
results in /-band and those in J, H and A'-bands in Lane et al] 
1 I2OOII) clearly indicate that Ba is brighter than Bb. This together 
with the I - J color favors a half subspectral class earlier for Ba 
than for Bb (see Fig.|6ll. 

The astrometric quality achieved with FastCam allows to lo- 
cate two new points on the GJ569Bb orbit around GJ569Ba. The 
orbits in the literature and the one derived including our points 
do not differ significantly and therefore the orbital parameters are 
in perfect agreement with those previously published, although 
our mass estimate of 0.116±0.007 M0 us ing the updated Hip - 
parcos parallax distance of 9.65±0.16 pc in Ivan LeeuwenI ( l2007h . 
Our mass estimate of the multiple GJ569B system is somewhat 
smaller, but within error bars of previously publishe d values ex- 
cept fo r the newly derived mass of 0.140^Q Qyg Mq in lDupuv et al.l 
ilOld) . Our WHT dat a point on June 20(39 falls within 1 sigma 
from the Dupuv et al ? mid) orbital solutions and our own or- 
bital solution although our NOT data point on July 2008 is more 
consis t ent with previous or bi tal solutions in 'Zapatero Osorio et al.l 
( I2OO4I) : ISimon et al.l ( l2006h : iKonopackv et aL (2010). For a sub- 
stantial refinement of the orbital parameters it would be necessary 
to sample the whole orbit with similar uncertainties as those derived 
from our observation at the WHT on June 2009. 

The data available on the GJ569B system is consistent with a 
primary of 0.081 ±0.010 Mq and a secondary of 0.059 ±0.007 Mq . 
If Ba were a binary system then the Bab component would have a 
mass < 0.018 M© . A simple qualitative analysis on the deviations 
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from a Keplerian orbit allows to place an upper limit to the product 
of mass and orbit semi-axis of this object. 
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